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THE ANISOTROPIC STRENGTH OF ASPHALT PAVING MIXTURES 


E. SHKLARSKY AND M. LIVNEH* 


Highway and Soil Engineering Laboratory, Technion-Israel Institute of Technology, 
Haifa 


ABSTRACT 


A laboratory study of compacted asphalt paving mixtures has shown that their strength 
is not isotropic. Preliminary test results are reported and discussed, and a modification of 
the bearing capacity formula is proposed. 


INTRODUCTION 


The bearing capacity of a compacted asphalt mixture depends primarily on its 
strength characteristics. Hitherto, bearing capacity predictions have been based on 
the assumption that this strength is isotropic, i.e. independent of orientation. This 
assumption, however, was not borne out in a laboratory study carried out by the 
authors. 

PRELIMINARY TEST RESULTS 


Two types of mixtures were prepared, one made with crushed limestone aggregate 
and the other with clean dune sand (see Figure 1). The bitumen used for both series 
was 80/100 penetration. The samples were compacted in 12 cm cubical moulds at 
an effort (energy per unit volume) equivalent to that of the Marshall method. 
Each type was studied by means of a simple compression test** in two series, with 
the load applied parallel and perpendicular to the direction of compaction respec- 
tively. The test results are plotted in Figure 2, and clearly show that a definite ratio 
m’, greater than 1, exists between the two respective strength values Gpar and Gperp: 


Eri SEs NG) (1) 
Operp 


The ratio was also found to increase with the density of the compacted asphalt 
up to a maximum value, after which it decreases with increasing density towards 
its asymptotic value of 1 (see Figure 3). It is believed that for the uncompacted 
asphalt the ratio equals 1. 


* The paper is based on the second author’s research work towards the D.Sc. degree under the 


supervision of the first author. 
** The asphalt concrete cubes were tested according to the Marshall procedure, i.e. after 


soaking in hot water (60°C) for half an hour; the sand-asphalt cubes were tested without 
preliminary soaking. The rate of compression was also according to the Marshall method, i.e. 2 


inch/min. 


Received November 18, 1960. 
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Figure 1 
Aggregate grading for bituminous concrete and sand-asphalt mixtures 


DISCUSSION 


The theory put forward by Casagrande and Carillo! on anisotropic materials is ex- 
tended here to include materials possessing both components of shear resistance, 
i.e. variable cohesion and constant internal friction. 


It is assumed that the compacted asphalt is anisotropic and has three strength 
parameters, namely y,—the constant angle of internal friction of the material for 
any given plane; Cpar—the cohesion coefficient parallel to the direction of com- 
paction; and Cperp —the cohesion coefficient perpendicular to the direction of com- 
paction. The value of the cohesion in a plane at an angle « with the principal hori- 
zontal plane (which is in turn at right angles to the direction of compaction) is, 
according to the above theory, 


G3 = Cocrp + (Cpar = Cperp) sin2a (2) 


If a compacted sample is subjected to the principal stresses operp (minor) and 
Gpar (major), the shear strength S under triaxial test conditions will be: 


Se Cperp aR (Cpar — Cperp) sin 2a = o,tang, (3) 


where &% is the angle of inclination of the critical plane, and oa, the stress normal 
to it. 
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Compressive stress at failure versus density (asphalt concrete cubes) 
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Figure 2b 


Compressive stress at failure versus density (sand-asphalt cubes) 
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nm’ versus density (points computed from Figure 2b) 


According to Mohr’s theory of rupture, the strength line given by eq. (3) is tangent 
to the critical stress circle. Hence (see App. I) the radius R of the latter is 


a (Cpa: + Operp tang,) tang, = 


+ /tan2 9,(Cpar + Operptanp,)? + (Cpa: + operp tangy,e)(Cperp + Operptan y,) (4)* 
The critical value of the angle a» is given by 


Hey fh atl Cperp + Operp tan Pe = 2Rtang, (5) 
9 Copar + Operp tang, 


For Cperp= Cpar, eq. (4) and (5) are reduced to the well-known equations of the 

critical stress circles for isotropic materials. 
As shown in App. II, eq. (4) can be rewritten with sufficient accuracy in the 
following form 
Ie Cpar tang, + J C2par tan2y, + Cpar Cperp + oper ee 


1 — sin 9, (6) 
The line S’, tangent to the critical circle, is thus a straight line of the form: 


S’=C+o,tan i. 
where C is given by id (7) 


C= (Coss tan Oats \/C2partan2p, + Cpar Coen) é oe a © 


* This expression can also be obtained by substituting y; =~3 =ye = const. in eq. (7) in Hank 
and McCarthy’s paper?. 
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Figure 4 
m versus n for g = 15°, 30°, 45° 


| S’ is in fact the strength line obtained in the triaxial compression test, C being the 
| apparent cohesion and g the angle of internal resistance. It is thus seen that C is 
| dependent on Cperp, Cpar and y,, and p equals p, provided eq. (6) is sufficiently 
| accurate. Eq. (8), in turn, can be rewritten in the following form: 


C= Cpoar (“sores | (tan Pe + [en Pe + 4 (9) 


COS , 


/ where 
Coar 
i 


(10) 


Coerp 


| As already stated, eq. (9) is valid if opar is the major principal stress; if it is the minor 
) principal stress, the apparent cohesion C is replaced by C’, which equals: 


pate. Cpar 1 “si sin Pe [Pa Bhar 
oe i ( COS @. ) (‘an Ye + /tan? yp, + "| (11) 
their ratio being 
| C n( tan De + taney. +1) 
| C tang, +./tan2o, +n 
| The ratio m is the counterpart of m’ [Eq. (1)] allowing for the difference between 


the simple compression test carried out on the asphalt cubes and the unconfined 
compression test. The existence of m’ also proves that of m, which is plotted against 
_nfor 9, = 15, 30, 45° resp. in Figure 4; the diagram explains the high values of 
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Bearing capacity [(a) and (b) — isotropic cohesion; (c) — anisotropic cohesion] 


m’' obtained for the asphalt concrete mixture compared with the low values obtained | 
for the sand-asphalt mixture: fora given n,m’ increases with increasing y,, and the 
latter is in fact higher for asphalt concrete mixtures. 


MODIFICATION OF THE BEARING CAPACITY EQUATIONS 


- McLeod’ suggested the following equation for the bearing capacity of asphalt 


paving mixtures: 
a 4C 1 + sin 
(pe eS ee) 


where p, is the ultimate pressure and C the apparent cohesion*. For the assumptions 
underlying this equation, see Figures 5a and 5b. 


* As the anisotropy of the material is not taken into consideration in McLeod’s equation, it is 
obvious that the apparent cohesion in question is the conventional value obtained in a triaxial test, 
with opar the major principal stress. 
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Figure 6 
Pm/pn versus m for py = 15°, 30°, 45° 


However, taking the anisotropic cohesion into consideration, eq. (13) becomes 
(see Figure Sc): 


nse 4C 1+sing_ ,l+sing : 
Pai isp ie ( m +1—sing) (lg) 
or 
1+ sing ; 
De px 4(—_— + 1 - sing) (15) 
The ratio = is plotted against m in Figure 6. For m=5 and 9 = 30°, a = 0.4. 


This result 1s of importance, since due to the anisotrop c cohesion the safety factoe 
for the allowable bearing pressure is considerably lower acording to eq. (13). Boths 
eq. (15) and eg. (13) have to be modified further to allow for factors not taken into 
account. In fact, eq. (15) is only an attempt to show the importance of considering 
the anisotropic strengths in the case of asphalt paving mixtutes. 


APPENDIX I 
According to Mohr’s theory, rupture occurs when the difference between strength 
(S) and shearing stress (t) is a minimum and zero. Under triaxial test conditions 
with a given radius R of the stress circle, S and t are functions of «, the angle of 
inclination of the plane under consideration. The rupture conditions are thus 


dS dt 
da dx 
S—t=0 (b) 


(a) 


and 
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Bearing in mind that t = Rsin 2a, eq. (a) becomes 


ae = 2Rcos a (c) 
da 
and eq,_(b): S = Rsin2« (d) 


Substituting R from eq. (d) in eq. (c), we obtain: 


dS 2S cos « aS 
‘da ce ada sin2Qe  tan2« 


(e) 
and substituting S from eq. (3), eq. (e) becomes 
2 cot Yeo [Cperp + Operp tan Oo Fr (Cpar = Cperp) sin2a% + 2Rcos2a, tan Ve| oS 


= (Cpar — Cperp) Sin 2%) — Rtan py, sin 2a. (f) 


Rearranging eq. (f) we obtain the critical value of the angle of inclination [eq. (5)], 
and finally, substituting eq. (5) in eq. (d) the radius of the critical stress circle [eq. (4)]. 


APPENDIX II 


The expression for R (qe. 4) can be rewritten as 


aa 7. . Btan Pe\ 2 
R -|( te Btang,) tan gt |(0 <7 Sone nF oa | Coat (a) 
where: 
ig Operp par 
p=, Jiang, +4; B= Pete GS Si rand 
n Cperp Cperp 


o, if 
A= tang,| 1+ 2tan2 —— Ye 2 nats 
of an. Petz COS P, tan v5) 


Neglecting the term BA, eq. (a) becomes: 


Btan, 


cos aa Car (b) 


Ree la + Btan,)tang, + D + 
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Percentage error E 


Figure 7 


Percentage error E versus B for p = 20°, 30°, 45° at n = 10 


which is identical with eq. (6). The percentage error E involved in omitting the term 
BA is 


B , B : 


eR. x 100 (c) 


E= 


Btan %, 
COS P, 


2 
(1+ Btang,)tang, +, |(D + ) + BA 


showing that E increases with n. 


E is plotted against B in Figure 7, for n = 10 (a reasonably high value) and prac- 
tical values of y,. The diagram shows that the maximum possible error under these 
circumstances is about 5% for m, = 20° and about 1.5% for p, = 30°. According 
to Nijboer*, the accuracy of the triaxial compression test is 30% for m and 10% 
for C. Substituting these values in eq. (a) we obtain 

1 
Eo 0.1(1 + Btan ¢,) (d) 
where E’ is the percentage error involved in predicting R in the triaxial compression 
test for isotropic materials; comparing it with the error plotted in Figure 7, it can 
be conclusively stated that the accuracy of eq. (b) is sufficient. 
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THE KINETICS OF THE ISOTHERMAL PEARLITE REACTION 
IN A 2% SILICON-MANGANESE STEEL 


S. NIEDZWIEDZ AND A. TAUB* 
Israel Institute of Metals, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


The transformaticn characteristics of a silicon steel were investigated in the pearlitic zone. 
The main cbject of the study was identification of the relevant area in the time-temperature 
transformation diagram and determination of the isothermal reaction. The activation 
energy of the process was calculated. 

The morphology and mechanical properties of the pearlites formed at different tempe- 
ratures were also observed. 


The pearlitic reaction of a high-silicon steel of the composition : C-0.63 %, Mn-0.99 %, 
Si-1.88%, Cr—0.41 %, S—0.006 %, P-0.037°%, was studied. The beginning and end 
of the transformation were determined under isothermal conditions for different 
| temperatures, whereby the pearlitic zone of the time-temperature transformation 
(TTT) diagram of this steel was identified. The transformed fractions were measured 
as a function of time under the microscope and the results used in the analysis of 
the kinetics of the reaction. The morphology and mechanical properties of the de- 
composed product were also observed. 


The procedure employed is described schematically in Figure 1. Small specimens 
are auStenitised (step a) and subsequently transferred (step b) to a constant-tempera- 
ture bath of liquid metal at temperature f, (between the A, and “knee” temperatures) 
at which they are kept for progressively increasing periods of time, resulting in 
' increasing amounts of pearlite in a decreasing amount of austenitic matrix at the 
given temperature. After the corresponding period, each group of three specimens 
is transferred very rapidly (step d) to a water quenching bath at room temperature, 
_ where the remaining austenite is transformed into martensite. The specimens thus 
| contain both pearlite and untempered martensite, differentiated under the microscope 
| through the selective action of the etching medium (nital or picral) which reveals 
them as dark- and light-etch areas respectively. 


i In the present experiments, specimens 8 mm in diameter and 2 mm thick were 
heated in a molten-tin bath, in three parallel series, at successive temperatures 
- from 700°C to 625°C in 25°C intervals. Each specimen was placed in the heating 
| chamber of a furnace for five minutes! at the austenitising temperature (830°C) 


* The Arturo Gruenebaum Chair of Mining and Metallurgy. 


| Received November 27, 1960 
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TEMPERATURE 


— 
aod 


TIME 


Figure | 
Schematic descripton of heat treatment 


and then transferred to the bath for the predetermined period of time and finally 
to the water bath held at room temperature. The treated specimens were subjected 
to microscopic examination in which the relative proportions of light and dark 
areas in the microstructure (as shown in Figures 2 and 3) were determined, using 


sa 


Figure 2 
Magnification: x 500 Biche Ricraly4ay, 
Specimen kept isothermally at 700°C for 40 seconds after austenitisation and subsequently quenched 
to room temperature, Dark areas 15 % pearlite: light areas 85 9 martensite 
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Figure 3 


Magnification: 500 Etch: Nital 2% 


Specimen kept isothermally at 650°C for 25 seconds after austenitisation and subsequently quenched 
tc room temperature. Dark areas 85°% pearlite: light areas 15°% martensite 


the electric semi-automatic counter developed by Hulburt,3 by Rosiwal’s method 
of lineal analysis?, based on the measurement of line intercepts across the different 
areas in a number of random directions (four in the present case, making a total 
of twelve measurements per series). * 


The mean, standard deviation, and standard error of the mean were calculated 
from the readings obtained for each series. The four traverses made gave a 95% 
probability of the error of the mean not exceeding + 1%. The average for each 
time-temperature combination was plotted in the diagram shown in Figure 4, giving 
the pearlite content (in percent) as a function of time, with temperature as a parameter. 
These curves are the actual reaction curves for the process. 


The reaction curves of Fisure 2 were used for the construction of the pearlitic zone 
of the TTT diagram of the steel in question, shown in. Figure 5. The graphical method 
employed is illustrated in Figure 6. In view of the difficulty in identifying the exact 
beginning and end of the reaction, the 1 % and 99 % transformation levels respectively 


were used instead. 


* Tn order to eliminate the operative error factor, the specimens were tested completely at random 
instead of in order of treatment. 
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Figure 4 
Experimental isothermal reaction curves, giving the pearlite content in percent as a function of time 
at various temperatures 
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Figure 5 
Method employed for construction of the TTT curve from the isothermal reaction curves 
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The rearlitic zone of the time-temperature transformation curve 
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} The sigmoid reaction curves of Figure 4 served as a basis for the determination 
| of the kinetics of the reaction. The curves are represented by the equation: 


} y=i1-e™ (1) 


where: y—fraction of transformed austenite; 
¢t — reaction time; 
m,n — constants for a given curve. 


—————————— 


Differentiating eq. (1) with respect to time, we obtain 


| dy pare n—1 j—mt" (2) 


| The maximum point of the derivative is the point of inflection of the sigmoid 
} curves. Differentiating eq. (2) and equating to zero, we obtain 


Pe n me Ac 
C2) ye [—nmt?""' + (n-1)t""7] = 0 


| dt? 

i whence nmt” = n-1 

| ‘eat (3) 
i substituting in eq. (1) ‘ots 
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The calculated values of n for the temperatures used in the present experiments 
are given in Table I below: 


TABLE I. 


Constants n,m for isothermal reactions 


Sopot 5 95/5 re eats limit m( x 10-8) 
700 4.43 +0.4 1.122 
675 4.64 +0.7 17.740 
650 ey) +0.4 4372 
625 3.16 +0.8 3457 


The table shows that around the “‘knee”’ of the TTT curve n= 3, for which y=1— et= 
= 0.45, while near the critical Ar, point n = 4, for which y = 1—e74 = 0.51. 


These results have a definite physical significance. Up to the point of inflection 
the reaction is accelerated; beyond it, it begins to decelerate due to the effect 
of impingement, while the point of inflection itself is characterised by an intermediate 
maximum rate of transformation. 


It can be seen that the points of inflection of the transformation curves occur at 
lower fractions of the total transformation as the transformation temperature ap- 
proaches that corresponding to the “knee” of the TTT curve. The constant n can 
thus be used as an indicator in comparing the rates of pearlitic transformation for 
different isothermal reactions. The variation of n is due to the formation morpho- 
logy of the pearlites which is due in turn to the rates of nucleation and growth. 


Eq. (2) can be used for calculating the activation energy of the process. In fact, 
it was shown by Zener‘ that the product nm provides for the temperature depen- 
dence of the reaction velocity, embodying as it does the rates of nucleation and 
growth (both of them temperature-dependent) conforming to the Boltzman distri- 
bution law, 


nm = Ce 22/RT (4) 


(Q, — activation energy, C — constant). Using eq. (1) in logarithmic form for a 
given value of y, we obtain 


1 
—In(i —y)=In a =mt"=b 


a 

: b 
where bis constant for a given case, or t = J —— 
m 
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Substituting m from eq. (4), we obtain 


ale "I nb 20,/RT 
Cc 


or t= C’e&/RT (5) 
where Q, is defined as Q,/n; or, again in logarithmic form, 


Int = —+1n C’ (6) 


If In t were plotted against 1/7, one would expect a straight line for Q,. However, 
in the pearlitic zone a curve with a negative slope would be obtained, since the 
graph would have to resemble a TTT curves. Hence in this case Q, can only be 
calculated from 


te = Cl eit TRVEE (7) 
where H is temperature-dependent. For the pearlitic-bainitic transition zone, 


Q, = 72,000 cal/mol; solving eq. (7) graphically for 700°C as shown in Figure 7, 
we find Q, = 77,000 cal/mol or Q, = nQ, = 308,000 cal/mol. 


inc! Q/RT 
Q\/RT Int, ] Inc? tnt 
dae ta ———— 
tines 


cat T 
PEARLITIC ZONE TRANSITION ZONE 
QO; H O;, 
Sn Cas ote rsa Int =InC’ + — 
Int = In C’ + Rr te RT n in Ete RT 
CONCLUSION: = tan w 
Figure 7 


Grarhical solution of eq. (7) for 


A peculiar feature of the steel under consideration, as observed under the micro- 
scope, is that a lamellar structure is also formed below the “knee” of the TTT curve 
at 625°C. This shows that the pearlite reaction reaches a maximum rate at the “knee”’. 
while lamellar structures are formed both above and below this point. Unfortu- 
nately, the methods used in the present study were not sufficiently sensitive to esta- 
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Figure 8 
Mechanical properties of the pearlitic structures as a function of the transformation temperature 


blish whether the rates of nucleation and growth also pass through a maximum 
and further research on this point would be justified. 


Below 625°C a transitory structure is formed, apparently a mixture of lamellar 
pearlite and feathery bainite. 


The steel under consideration was found to be nearly eutectoid. Only traces of 
hypo-eutectoid ferrite could be detected (precipitated at the former austenitic grain 
boundaries) and their influence could only be negligible. 


Finally, the transformed specimens were studied with a view to determining the 
effect of variations in the inter-lamellar spacing of the pearlite upon the mechanical 
properties of the steel. Hardness tests were carried out on the Rockwell C scale (limit 
of accuracy +1 point) and tensile tests ——on a Hounsfield tensometer, using 
specimens 3.57 mm in diameter and 12.6 mm in length, i.e., with an initial cross- 
section of 10 sq. mm and length-to-diameter ratio of 4. The average elongation 
for the entire pearlitic zone was found to be 15%, witha slight increase at the “‘knee’’, 
and the Charpy impact tests showed an average value of 0.4 kg-m. Plastic deformation 
prior to “necking” was relatively small, which indicates a low value for the coefficient 
of strain hardening. The relation between hardness, U.T.S. and transformation 
temperature is shown in Figure 8. 


The machinability of all the pearlitic structures of the steel was found to be low 
(in spite of the low coefficient of strain hardening); hence it is not recommended 
for operations involving extensive machining. 


Comparison of the hardness values reported here with those obtained on tempering 
martensitic structures of the same steel’ reveals that isothermally formed pearlite 
is equal in hardness to the decomposition product obtained by 30 minutes’ tempering 
at similar temperatures. 
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THEORETICAL ANALYSIS OF THE SPLITTING TEST FOR ASPHALT 
SPECIMENS 


E. SHKLARSKY AND M. LIVNEH* 


Highway and Soil Engineering Laboratory, Technion-Israel Institute of Technology, 
Haifa 


ABSTRACT 


The splitting test used in determining the strength of concrete specimens is applied to | 
asphalt paving mixtures. A theory is proposed for the failure mechanism, using a stress _ 
distribution obeying Hooke’s law for the elastic stage, and an anisotropic strength curve | 
assuming a constant angle of internal friction and a variable cohesion coefficient—for the | 
plastic stage. 


Satisfactory agreement was obtained between the proposed theory and laboratory | 
experiments. The importance of these findings lies in the possibility of utilising the theore- | 
tical equations in determining the relationships between the failure stress and the height | 
of the specimen, as well as the anisotropic cohesion from perpendicular and parallel split- | 
ting tests. | 

| 


1. INTRODUCTION 


The ideal procedure for determining the strength of asphalt specimens is one which 
would permit laboratory design of the product on the one hand and field control © 
on the other, using an analytical method based on strength properties governed 

by the quality of the product but independent of the type of test. This requirement 

is only partially satisfied in the testing procedures in use at present. For example, | 
the Marshall test! permits both laboratory design and field control, but without | 
any reference to the strength properties described above. Conversely, the triaxial | 
shear test!, although based on the angle of internal friction and cohesion, does | 
not permit field control.** 


* The paper is based on part on the second author’s research work for the D.Sc. degree under 
the supervision of the first author. 

** Field control of asphalt paving mixtures by the triaxial test technique is impracticable owing _ 
to the requirement limiting the height-diameter ratio of the specimen to 1.5—2. The thickness of | 
the asphalt layer being about 5 cm, the corresponding diameter would be 2.5-3 cm; such a specimen | 
would necessarily be unrepresentative considering the grain size of the asphalt( <2 cm), the minimum 
required ratio between diameter and maximum grain being 4 : 1. bg 


Received May 9, 1961. 


202 Bull. Res. Counc. of Israel, Vol. 9C, 1961 


Vol. 9C, 1961 E. SHKLARSKY AND M. LIVNEH 203 


These circumstances have led the authors to adopt a different method constituting 
a far closer approach to the ideal testing procedure. The splitting test*, also known 
as the “‘Brazilian test’? and commonly used in determining the tensile strength 
of concrete, consists in applying a load by means of two parallel sticks to a con- 


crete cylinder along opposite generatrices, or to a prism along the centre lines of 
opposite faces’. 


The splitting test has also been adopted for soils*~°; its application to asphalt, 
to the best of the authors’ knowledge, is proposed for the first time. The recommen- 


dation of a prismatic specimen is based on considerations of field control require- 


ments, the undisturbed specimens having to be prismatic in order to lend themselves 
both to horizontal and to vertical loading. 


2. THE ELASTIC LIMIT STRESS IN THE SPLITTING TEST 


The theory underlying the splitting test for concrete is based on a single property, 
namely tensile strength. This approach is impracticable for asphalt in view of its 
failure to agree with test results (see Section 5). In these circumstances, the proposed 
analysis is based on the characteristic strength properties of asphalt, namely the 
angle of internal friction, w, and cohesion coefficient, C. 

The elastic stress distribution along the diameter of a cylindrical specimen 
perpendicular to the loading stick is given (assuming the material to be weightless), 
according to Peltier’, by eqs. (1a)—(1c) (see Figure 1). 


2pa | R2 ie 

= Pee (5) 22 = la 
pe le (6 — $sin 26) | (la) 

2pa | R? “lar 2R 
eee ee —1 1b 
¥ oR | 6 + $sin28) + (1b) 
Wee — @) (1c) 

a 
Looaing ae ve) 


orR 
2ap 


2a 


Figure 1 le 
Elastic stress distribution in z-direction in splitting test 


* For a literature survey on the splitting test for concrete, see Ref. (2). 
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where R is the specimen radius, 2a — the width of the loading sticks, p — the com- 
pressive stress uniformly distributed across each stick, o, and o, — axial stresses 
in the x- and z-directions respectively, t,, — shearing stress in the xz-directions, 
and the angle 6 is defined as follows: 


tand = A (id) 


In accordance with the above equations a, is compressive* and o, compressive; 
at the edges of the specimen and tensile within it. 
Assuming that the strength of the asphalt obeys Coulomb’s law, namely that: 


S=C+4,tang (2}) 


| 
where S is the shearing strength, C — the cohesion coefficient, ~ — the angle of 
internal friction and o,—the stress normal to the plane of failure, then the following: 
condition has to be satisfied at every point of the plane of failure: | 


6,-0 C cos o, sin | 
( 1 3)f fal Y 3 Pp (3) 
| 


2 1—sing 1-—sing 


In this equation (o, — 03), denotes the difference of the principal stresses at failure 
(o, being the major and o, the minor principal stress). In practice, o, in eq. (1) 
represents o, and o, represents o,, so that substitution of their difference in} 
eq. (3), giving eq. (4), permits determination of the stress py corresponding to the: 
limit load in the elastic range:** 


p x 2a a pee: R _ Cocos 
2 read, 


mR az ~ L-sing 
Dix 2a RR Wee sing 
Siren E (6 — $sin 26) — 1 Trai’ (4) 
Eq. (4) can also be rewritten in the following form, designating EIN C8? by F: 
pa 
= (cas 4 (sin26 — 26sing) + —{———_x21—sing)+2sing (5) 
27 tand + 1 | 


* In this paper, the plus sign denotes compression and the minus sign tension. 
** Tt should be borne in mind that if this limit is exceeded, stress distribution will be affected andt 
eq. (1) will no longer be valid for calculating the failure load. 
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The value 69, for which F is maximum, determines the coordinates of the initial 
point of the failure mechanism, namely for 


dF 
rene (6) 
we obtain 
sin 269 — 269Si 1—si R 
- eer ge ee de ae ( 3 Ff ra 
; (2—tand)+1)?acos”5, athe A Ve 
a 
whose approximate solution* is 
tan dp = tania af) ae garg 
RARE E (8) 
WN 


where 


1 2 
No = tan? (1/4 — 9/2) = tan (2/4 ae y/2) 
It can be shown that for R/a = 00, the solution according to eq. (8) is quite exact. 
Finally, substituting eq. (8) in eq. (5), we obtain 


R — a 
CxR cos ae naa Ree Lteteey 
See = 2 ; s~sing arctan ( | + 
E ing Ae, JN 
1- 2 /Ng | No +(1-4VNe) JNe 
Dees: 
(i — sing) i 
2—-— JNo 


It is of interest to consider the case R/a = © (which is identical in practice with 
that of a semi-infinite elastic body, subjected to an infinite strip of load) as solved 
by Boussinesq. Substituting R/a = oo in eq. (9), we obtain 


Cn cos p 
DRI = IS Ea ae, (10) 


cosy — FS — p)sing 


* For example, for R/a =5 and 9 = 30°, the exact solution of do by eq. (7) gives 23°48’, while 
approximation according to eq. (8) gives 20°42’, the error involved in this approximation being 
(3.94-3.90)/3.90 =1% for F and (3.44-3.27)/3.27=5 % for Zo (the coordinate of the point of failure). 
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_A similar expression can obviously be arrived at by substituting the difference 
(co, —o,) given in Boussinesq’s formulae® in eq. (3) and proceeding as described 
above. 

The next step is to see to what extent the above analysis is applicable to the pris- 
matic specimens on which the proposed method is based. For 2a =0, solutions 
exist for: 


(a) A prism of square cross-section (solution by means of the relaxation method). 
(b) A prism with an infinite rectangular horizontal cross-section®*. 


Comparison of these solutions with one for a cylinder with 2a = 0 leads to the. 
conclusion that stress magnitudes and distributions are almost identical in all three | 
cases. | 

For example, the respective stress for each section, at the mid-point of the axis | 
joining the two points of application of the load, is given by 


cam 


C¢= 


A ry 
— 
— 
— 
See” 


2 


where P is the load per unit length, R — the radius or half-height of the specimen, 
and k a constant factor given in the following table: 


k for different sections 


Section 


Stress 


Infinite Square Circular 
rectangular 


— 0.484 
+ 1.824 
+ 2.308 


The numerical values in the table show that the elastic solution for the circular | 
section is very close to the average of the square and the infinite-rectangular. Hence | 
the conclusions arrived at above are also valid for the prismatic specimen. These | 
conclusions are: | 


* This solution is based on the superposition of two identical simply-supported beams of in- 
finite span, each subjected to the same single force at midspan acting in opposite directions, and | 
is obtained by means of integration of series. 


Vol. 9C, 1961 E. SHKLARSKY AND M. LIVNEH 207 


(a) The relative coordinates of the initial point of failure are independent of R/a, 
since by eq. (8), we have 
a 1 
Re == 


- 


“"(b) The elastic limit stress increases with R/a [eq. (9)] tending to an asymptotic 
value for R/a = o (eq. 10). 


3. THE FAILURE STRESS IN THE SPLITTING TEST 


The elastic limit stress pp is not identical with the failure stress p,. The former 
represents the initial point of the development of the failure mechanism and the 
latter — its end-point. 

From conclusion (a) in the preceding section, it follows that the failure mecha- 
nism in its early stages is the same for prismatic specimens, irrespective of the h/a 
ratio (h being the half-height of the specimen); if the initial point coincides with 
the centre of the specimen, the possible mechanism is that of unconfined compression'® 

From conclusion (b), it follows that the failure stress increases with the h/a ratio 
and tends to anasymptotic value for h/a = «.This asymptotic value is in practice 
the failure stress according to Prandtl’s theory®. Hence the failure mechanism in 
the splitting test has to correspond to either alternative according to the h/a range 
in question. Such a mechanism is described in Figure 2 (also giving the actual failure 
patterns obtained on asphalt specimens); for 0 =0 it corresponds to the uncon- 
fined compression alternative, and for 0 = 2/2 — to Prandtl’s theory. The area 
bounded by points 1, 2, 3, 4, 5 comprises the zone of Prandtl’s mechanism, lines 
5-6 and 2-6 being assumed as circular in their upper part and straight in the lower 
part,* forming an angle of (x/4—q/2) with the vertical axis of symmetry, which is 
the locus of the principal stresses. The transition point between the two parts is 
determined by the consideration that for 0=72/2, h/a = oo. This hypothesis is 
comparable to the failure mechanism proposed, on the basis of exact analysis, for 
a material with a curvilinear strength diagram’' (see Figure 3). 

The geometrical relationships for the proposed failure mechanism are 


re =roe"*?; x =rosin(x/4 + y/2); a =rocos(n/4 + y/2). (12) 


i) 
(1+ e @ sin 


~ sin (a/4 — 6/2)sin(x/4 — apices 7/2 + sing/2./No] + VNo (13) 


h/a 


* The combination of circular arc and straight line was assumed for simplicity. 


Figure 2a 


Theoretical failure mechanism in splitting te| 


Figure 2b | 
Failure patterns in asphalt specimens (2a = |! 
(1) 2h = 6cm 2) 2h = 12c 


Figure 2a 


Figure 2b(1) Figure 2b(2) 
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(a) (b) 
Figure 3(a) Figure 3(b) 
Critical curve of shear envelope— Failure mechanism for shear envelope in (a) 


curvilinear strength diagram (Ref. 11) 
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Figure 4 
Forces acting on segment of failure mechanism 


a 
< 
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The failure stress py can be determined by substituting the forces acting on a seg- 
ment of the proposed mechanism (see Figure 4). The sum of moments with respect 
to point 3 is 


2 2 t) ei 
Fa, zs pra = Epre i: Creat do (14) 
2 2 2 0 
where 
ppg me Os (15) 
N —— 
r ANG 
and 
E,=1,Ccos@ (16)* 


Solving the integral in eq. (14) and substituting the appropriate values from eq. 
(12), (15) and (16) we have: 


pr = Ccotp[e”?"?N,— 1] (17) 


Eqs. (13) and (17) give the relationship between p, and h/a with the angle @ as 
parameter. It can be seen that: 


(a) For 6 = 0, by eq. (13) h/a =./No and p,, according to eq. (17), is the un- 
confined compression. 


(b) For 0 = 2/2, h/a = «© and p, is the Prandtl compression. 


4. THE EFFECT OF THE ANIS@¥ROPIC COHESION ON THE FAILURE STRESS 


The strength of the asphalt being anisotropic, the effect of this anisotropy on the 
failure stress in the splitting test should be studied. The strength diagram for asphalt 
Ss ales 
is 


S = Coerp + (Coar — Cperp) Sin20%q + (6 par — Tperp) COS2% tan p (18) 


* The force Eg can only be determined provided 03, acting in the plane of failure 5-6 or 2-4, 
is known. Since this plane of failure does not forman angle of (/4-¢/2) with the horizontal boun- 
dary line (for which it is known that o3 is normal and zero), it must be assumed that the zone bet- 
ween the boundary and the plane forming an angle of (/4—9/2) with 5-6 or 2-4 is stress-free; hence 
63 for the above plane of failure is zero. 
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6 par 


6 perp = . f— — —>Coerp Overp 


Direction of compaction 


6 par 


Figure 5 
Scheme of principal stresses and cohesion relative to direction of compaction 


where (see Figure 5) 
Cyerp = Cohesion coefficient perpendicular to direction of compaction 


= cohesion coefficient parallel to direction of compaction 


par 


g = angle of internal friction (constant for each direction) 

Opar = Major principal stress, parallel to direction of compaction 

Oyerp = Minor principal stress, perpendicular to direction of compaction 

% = angle of inclination of plane of failure with respect to minor principal 


stress, given by 


ae Coar + Fperp tan Y (19) 


The strength diagram according to eq. (18) cuts the Mohr circles, while the tangent 
is given with sufficient accuracy by’? 


S’ =C,, + (0, — 03) cos? (2/4 + p/2) tan p (20) 
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where C,, is the apparent cohesion related to Crop, Cpar and p and the angle between _ 
the major (o,) or minor (03) principal stress and the direction of compaction. For 
the case where o,,, is the major stress we have C = C,,, and for that of the minor — 
stress, C,,= C/m; for intermediate cases the following interpolation can be used 


Cy»=C+ (= — c sin’ B (21) | 


where f is the angle between the major principal stress and the direction of com- 
paction. 

In the case of anisotropic cohesion, the failure stress will be determined on the | 
assumption of a strength diagram according to eq. (20). The principle underlying | 
this assumption is based on the common soil engineering practice (known as the | 
‘‘m = Oanalysis’’!>; see also Appendix) of approximating the actual stre~ eth diagram 
and failure mechanism by their apparent counterparts. 

The apparent mechanism corresponding to the apparent strength diagram [eq. 
(20) ]is one in which the angle between the failure lines and the principal stresses is 
m/4 + m/2. It is similar to that described in Figures 2 and 4, except that the cor- 
responding forces in this case are 


Es | xe 2 | x (22) 
Q 


E, = E ~ (<-c) sin | rg COS P (23) 


From the sum of moments with respect to point 3 we have 


E 4x Ppar a? E ro an 
5 + Saag a + i [C + (C/m — C)sin’]r3 e 7° '"? do (24) 
Solving the integral in eq. (24) and substituting the appropriate values from eqs. 
(12), (22) and (23), we obtain 


ar an 1 
Peer = [e??'? N, — 1] cot E pia 6 ce sing)sin*o | + 
(1 1 
e errant + sing) cos 0(e*°'"? cos @ — 2e79'"? sin@ —cos0) (25) 
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Eqs. (25) and (13) give the relationship between Pror and h/a with the angle 6 
as parameter (see Figure 6). For 0 = 0, eq. (25) is identical with eq. (17) which is 
its counterpart for an isotropic material, while for 0 = 7/2, eq. (25) is equivalent 
to eq. (17) multiplied by 4[1— sing + (1+ sing)/m] which is actually the aniso- 
tropy factor in compression according to Prandtl*. 


Eq. (25) permits calculation of the ratio Pyar/Pperp» given graphically in Figure 7 
for m = 3. The figure shows that this ratio is smaller than m. 


5. EXPERIMENTAL 


On the strength of the above, a number of theoretical relationships can be derived 
and verified in laboratory experiments These are: 


(a) The relationship between the respective failure stresses under parallel loading 
(p{° ) and the unconfined compressive stress (also parallel, o,), for h/a = 6; results 
given in Figure 8 for m = 1 and 3 respectively. 


(b) The relationship between the respective failure stresses under parallel loading 
for h/a = 6 nd h/a = 12 (pS and p&,?); results also given in Figure 8 for m=1 
and 3 respectively. 


(c) The relationship between the respective failure stresses under parallel and per- 


pendicular loading for h/a = 12 (p{i?? and p{}”) ); results given in Figure 9. 


Specimens were prepared from crushed limestone (for grading see following table) 
bitumen (6% by weight), subjected to static compaction at 140°C to give different 
densities and tested after half an hour’s soaking in hot water (60°C). 


Grading of aggregates 


AS.T.M. sieve 3/4” 1/2” 3/8” 4% 10% 4o* go* 2007 


% passing 100 80 68 sy 38 17 11 6.5 


Tests were carried out, at a steady rate of loading of 2 inch/min, as follows: 


(a) Splitting test, using a 6 cm specimen and Icm loading stick (h/a = 6), paral- 


lel loading (p\%?). 


* It is interesting to note that this anisotropy factor is identical with the one obtained when 
calculating the bearing capacity according to McLeod !2. 
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Failure stress/cohesion ratio vs. h/a (parallel loading) for m = 1 (solid line) and m = 3 (dashed line) 
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Failure stress ratio (parallel and perpendicular loading resp.) vs. 


h/a for m 
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(a) Failure stress/unconfined compressive stress ratio vs. angle of internal friction for h/a = 6 
(solid line) (parallel loading, m=1 and 3 resp.) 


(b) Failure stress ratio (h/a = 6 and 12 resp.) vs. angle of internal friction (dashed line) (parallel 
loading, m = 1 and 3 resp.) 
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Figure 9 


Failure stress ratio (parallel and perpendicular loading resp.) vs. apparent anisotropic cohesion, 
for h/a = 12 
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(b) Splitting test, using a 12cm specimen and the same loading stick as before 


(h/a = 12), lo ding perpendicular and parallel (p{t7) and p,ur 


(c) Unconfined compression test, using a specimen 4” in diameter and 6” high, 
parallel loading (¢,). 


(d) Compression test, using a 12cm cube, perpendicular and parallel loading 
(Greens Oper): 


Results are summarised in Figures 10, 11, and 12, giving failure stress against density. 
The ratios calculated as described at the beginning of this section are given in 
Figure 13. 

Comparison of experimental and theoretical data (Figures 8 and 9) showed satis- 
factory agreement. Values of y used for comparison should be above 40°. (It 
should also be borne in mind that the ratio m’ = Gy4,/Gperp iS equivalent to m in 
order of magnitude only). 

Ratios are summarised in the following table, showing that except for the low 
density of 2220 kg/m the deviation is of a reasonable order of magnitude, its maxi- 
mum being 20% for the p97/ p\2”) ratio and the 2240 kg/m} density. (For the 
2220 kg/m? density the maximum deviation is 27%, for the same ratio). 
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iF PP Io. EUS) WPL IANS PRT IRA MOXIE POLS. || Bil I Ailey POO | Dsee | 
Cas Ve 1.70 1.86 1.92 1.65 


3. pO?p?)) 1.76 | 1.29 | 1.38 | 1.65 | 1.33 | 1.42 | 1.54 | 1.34 | 1.44 | 1.44 | 1.24 | 1.32 


4. pi) ipO | 0.60 | 0.64 | 0.58 | 0.60 | 0.65 | 0.59 | 0.61 | 0.66 | 0.60 | 0.64 | 0.64 | 0.58 


6. CONCLUSIONS 


The theory of the failure mechanism presented in this paper is especially useful 
for analysing results obtained in splitting tests of field specimens of varying thick- 
ness. In addition, determination of the bearing capacity of the asphalt involves 
that of the m ratio, which is a function of the anisotropy of cohesion in field spe- 


cimens; but since only the p,4/Pperp Tatio lends itself to laboratory determination 
the relationship between m and the latter must also be ascertained. 
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Figure 10 
Failure stress vs. density in splitting test (h/a=6) and unconfined compression test, parallel loading 
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Figure 12 
Failure stress vs. density in compression test, parallel and perpendicular loading 
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Ratios obtained in splitting and compression tests vs. density 
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In view of the satisfactory agreement between theoretical and experimental re- 
sults, it can be concluded that: 


(a) The theoretical relationship between p,,, and h/a (Figure 6) can be used for 
calculating the equivalent failure stress for a specimen of standard h/a (recom- 
mended value — 6) on the basis of failure stresses determined on specimens of 
different h/a. 


(b) The theoretical relationship between the ratio P,o-/Pperp and h/a (Figure 7) 
can be used for calculating m on the basis of splitting tests under perpendicular 
and parallel loading. 


APPENDIX 


The failure mechanism for h/a = co in the case of anisotropic cohesion and zero 
internal friction angle is described in Figure 14a. 
The forces involved are (see Figure 14b): 


Cre Cet (C. = C.cysineg 
E, — (p- 2eCos Gras) x (a)* 


E,= 2./ Coar C perp X 


The geometrical relationships in this mechanism are 


X= atana,; b=a/cosda); d= 2a/tan24 (b)** 


r= a/sinag; .tandg =) Oya pl Cour 
The sum of moments with respect to point 0 gives 


2 to 
x a Re 2 
ae aban oie — = 0 Cc 
EF +4 }erS E, (5 +4) +2f Cds (c) 


* Obtained by substituting ¢ = 0 in eq. (18). The radius of Mohr’s circle described by varying 


@ is 4/ Coerp Cpar (See Figure 14c). 


** Obtained by substituting ¢ = 0 in eq. (19). 
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Figure 14 


(a) Failure mechanism in semi-infinite specimen 
(b) Forces in failure mechanism (a) 
(c) Anisotropic envelope for pg = 0 


Solving the integral in eq. (c) and substituting the appropriate values from (a) 
and (b) we have 


age 4(2C erp C mee ar = 2(C erp C a, er. +2(C er. +C an 
p= (CrarlCm | perp P. P Ge oe (Cocep + Cy (d) 
perp par) “ par 


Cc 
tan% = % = ee for Cpar!C perp > 0 (e) 
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hence the bearing capacity is 


14 Cie Cree = OC ia: aN, Peli eek ie 
(3C ap a Char) C var V Car Coerp 


C far eSr Do 
and for —"" = o, pS Ax Ce Crerp (f) 


perp 


The apparent strength diagram of the asphalt* is given by 


Cia) CoCr, == CONSt. (g) 
G=.0 
to = 1/4 


hence the bearing capacity calculated from it is 


pie 512 */ Cry Coe 


Comparing (h) with (f), it can be seen that the error involved in substituting the 


apparent strength line is only 10%, and this only for the extreme case of C,,,,/C 


perp 
oo. For smaller values the error will decrease, hence the use of this simplified 


method is justified in practice. 
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A SYMPOSIUM ON BASIC PROBLEMS OF SPACE TRAVEL 


Monday afternoon, 3.4.61 Chairman: D. ABIR 


Interplanetary matter and planetary atmospheres 
J. NEUMANN, Department of Meteorology, The Hebrew University of Jerusalem 


This paper reviews progress since the launching date of the first artificial Earth 
satellite. 

Interplanetary matter. The major discovery is the existence of the Van Allen 
radiation belts around the Earth extending outward to several Earth radii. There are 
some indications that, at least, one of the other planets (viz. Jupiter) possesses 
similar belts. 

The atmospheres of both the Sun and Earth (and, presumably, those of other 
planets as well) extend to greater distances than assumed heretofore. The Earth, 
e.g., is apparently within the Sun’s outer corona. Influx of matter fiom solar storms 
causes pronounced variations in upper atmospheric densities. Data have been gathered 
on meteotitic particles about the Earth and on the rate of accretion to Earth of 
such matter. Clouds of ionized particles, cairying a magnetic field, are likely to 
be encountered in interplanetary space. 

Planetary atmospheres. A considerable body of data has been obtained on the 
outer reaches of the Earth’s atmosphere, while special satellites have taken photo- 
graphs showing that cloud systems in the lower terrestrial atmosphere have a spiral 
structure. The latter give rise to new ideas on circulation conditions within these 
storms. Considerable work has been done on the subject of the tenuous lunar atmos- 
phere. 

Concerning the other planets, the most significant discovery is the definite identifica- 
tion of (small amounts of) water vapour above the visible surface of Venus. The disco- 
very was made from a balloon, without recourse to the more expensive rockets or 
satellites. Earlier ideas, suggesting a thick atmosphere over Venus, are giving way to 
hypotheses supporting a relatively thin atmosphere. Because or the high planetary 
albedo of Venus, its average atmospheric temperature is probably lower than that 
for Earth. Fundamental astronomic data in respect to Venus are still uncertain, 
but if it is true that the rotation period of the planet is something like three (terres- 
trial) weeks, then, it is probable that meridional circulations dominate the atmos- 
pheric motion. 

Little additional knowledge has been gained in respect to the other planets, but 
indications of a low-level plant life on Mars have strengthened. 
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224 BASIC PROBLEMS OF SPACE TRAVEL 


Reentry of space vehicle into planetary atmosphere 
J. Rom (Rasinowicz), Department of Aeronautical Engineering, Technion—Israel 
Institute of Technology, Haifa 

The flight conditions of a space vehicle reentering a planetary atmosphere at 
extremely high speeds (18,000-33,000 ft/sec) have to be known for the design of a 
recoverable manned or instrumented space vehicle. The extremely high energy 
that is dissipated in the deceleration of the vehicle is manifested in the heating of the 
vehicle end in thermodynamical changes in the state of the air. Some of these prob- 
lems and suggested solutions will be discussed. 


Terrestial radiation belts 

H. Kasua, Physics Department, Technion-Israel Institute of Technology, Haifa 
An account is given of the discovery of the Van Allen radiation belts and of the 

evidence obtained in flight of satellites and rockets, on the form of the belts and 

the nature of the particles trapped in them. The implications of the discovery on 

manned space flight is discussed. 


Basic problems of space travel 
A. YARON, Scientific Department, Ministry of Defense 

a. Theorem of conservation of momentum; rockets as means for obtaining 
high flight velocity; stage rockets and staging techniques. 

b. Rocket efficiency parameters; specific impulse; propellant characteristics; 
review of the state of the art. 

c. Design problems; search after construction materials for space travel. 

d. Problems of human space travel; weightlessness, resistance to acceleration. 
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